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ABSTRACT 

A new method is introduced for making X-ray mass determinations of spherical clusters of galaxies. 
Treating the distribution of gravitating matter as piecewise constant and the cluster atmosphere as 
piecewise isothermal, X-ray spectra of a hydrostatic atmosphere are determined up to a single overall 
normalizing factor. In contrast to more conventional approaches, this method relies on the minimum 
of assumptions, apart from the conditions of hydrostatic equilibrium and spherical symmetry. The 
(3JT). method has been implemented as an XSPEC mixing model called CLMASS, which was used to determine 

masses for a sample of nine relaxed X-ray clusters. Compared to conventional mass determinations, 
CLMASS provides weak constraints on values of M500, reflecting the quality of current X-ray data 
for cluster regions beyond r^g. At smaller radii, where there are high quality X-ray spectra inside 
and outside the radius of interest to constrain the mass, CLMASS gives confidence ranges for M2500 
that are only moderately less restrictive than those from more familiar mass determination methods. 
The CLMASS model provides some advantages over other methods and should prove useful for mass 
determinations in regions where there are high quality X-ray data. 

Subject headings: galaxies clusters: general — intergalactic medium — methods: data analysis - 
k«J ■ X-rays: galaxies: clusters 
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1. INTRODUCTION 



The need for accurate masses for galaxy clusters is motivated by more than their significance as fundamental 
parameters of the most massive, dynamically stable structures in the universe. Because gravity is dominant o n the 
£2 ' large st scales in the un iverse, masses of the larg est virialized structures are important cosmological probes (e.g. iVoitl 
120051: lAllen et~al . 2008; Samushia & Ratra 2003). For example, the growth of structure at high mass scales provides 
a sensitive means to measure the power spectrum of density fluctuations that emerged from the early universe and 
to deter mine the time depende nce of the rate of cosmic expansion, hence placing useful constraints on cosmological 
> ; models (jVikhlinin et al.ll2009b| ). 

The dominance of gravity in clusters also means th at there are severa l avenues for clus ters mass determina- 
0^ ■ tions, which can be classified broadly as dyna mical (e.g. iGirardi et alj|1998t iRines et al.ll2003D. gas d ynamical (e.g. 
; iPointecouteau et al.l2005t I Vikhlinin et al.ll2006[ ) and lensing fe.g. lSmith et al.l2005HLimousiri et al.l2007ft . No approach 
is free of weaknesses. Dynamical mass deter minations are limited by our inability to measure velocities transverse to 
OO ' the line-of-sight (e.g. iMamon fc Bouc 2009). Gas dynamical methods rely on the intracluster medium (ICM) being 
f^) ' hydrostat ic fcf. | Markevitch fc Vikhlininl l2007fl and t hev require accurate determinations of the total pressure of the 
f^*) | ICM (e.g. IChurazov et al.ll2008 | : | Mahdavi et al.ll2008[ ). Strong lensing is confined to the highest density regions of rich 
clusters (e.g. Newman et al.ll2009f ) . while weak lensing is affected by substantial statistical uncertainties for individual 



clust e rs, line-of-sight pro jections and the mass-sheet degeneracy (e.g. iMetzler et aHl2001t lHockstra 2003; Bradac fet al.l 
[200l iZhang et al.[|200l . Many of these issues can be addressed by judicious choice of targets, but the most reliable 
. results are likely to come from the application of multiple mass determination m ethods to cross check and augment 
j_j ■ one another (e.g. iMahdavi et alJl20"07t IChurazov et alJl2008t iNewman et~aT1l2009H . 

C3 1 This paper introduces a new approach to X-ray mass determinations. The method relies on the same basic as- 
sumptions as other X-ray mass determinations, i.e., that the intracluster gas is hydrostatic and spherically symmetric. 
However, it largely avoids the additional, model-dependent assumptions that are required by other approaches. The 
method is outlined in section [2] and applied to a sample of clusters with high quality Chandra data in section |3l 
Strengths a nd weaknesses o f the m ethod are related to the results in section |H 

Following iVikhlinin et al.l (|2006f ). all distances are computed assuming a flat ACDM cosmology, with J7m = 0.3 and 
a Hubble constant of 72 km s^ 1 Mpc^ 1 . For each cluster, the radii r$oo and r25oo enclose regions with densities of 500 
and 2500, respectively, times the critical density at the cluster redshift, while M500 and M2500 are the corresponding 
enclosed masses. Confidence ranges are determined at the 90% level except where stated otherwise. 

2. METHOD 

X-ray mass determinations rely on the ability of X-ray observations to determine the temperature and density of 
the hot ICM. In an ICM dominated by hot gas, the pressure is determined by its temperature and density. Unde- 
tected components of the ICM, i.e., nonthermal particles, magnetic fields and turbulence, can add to the effective 
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pressure, although simulations (e.g. iNagai et' al. 2007) and observations (e.g. iChurazov et al.|[2008T > suggest that these 
contributions are modest, of the order of 10%, in relaxed appearing clusters. Since the purpose of this paper is to 
introduce a new approach to X-ray mass determinations, such effects are not considered further here, but they need 
to be understood for precise mass determinations. 

If the gravitational potential is spherically symmetric and the ICM is in hydrostatic equilibrium, the distribution of 
the ICM is governed by the equation of hydrostatic equilibrium, 

dp GM(r) 

= -Pn( r )—^—> (1) 

where p(r) is the gas pressure, p g (r) is its density and M(r) is the total gravitating mass within the radius, r. If the 
pressure and density distributions of the ICM are known, M(r) can be determined from this equation. In practice, 
properties of the ICM (particularly its temperature) are only determi ned at discrete radii, so that the pressure derivativ e 
is constrained poorly, if at all, by the data. With few exceptions dNulsen fc Bohringerlll995t lArabadiis et al.l feOO^. 
mass determinations circumvent this difficulty by employing model-dependent assumptions to supplement equation 
(fT). Details vary widely, but typical approac hes assume analytic forms, with a modest numbe r of parameters, eithe r 
for the mass profile (e.g. IMahdavi et al.ll2l)l)~i. or for the temperature and density profiles (e.g. IVikhlinin et aLl feOOB). 
For example, using the latter approach, substituting the analytic form for the pressure and density profiles in equation 
((T|) gives an analytic expression for (1/ p & )dp/dr, with parameters that can be determined by fitting the analytic 
temperature and density profiles to the discrete tempe ratures and densities obtained from X-ray spectra. 

Also with few exceptions (e.g. IMahdavi et al. I [20 071. several steps are required to obtain a cluster mass, so that 
a confidence range for the mass is determined indirectly from fitting the X-ray spectra. For example, data may be 
"deprojected" to obtain gas temperatures and densities (with strongly correlated errors). Deprojected temperatures 
and densities are then fitted with analytic models to determine their parameters, which are used to obtain the mass, 
as described above. Doing mass determinations in several steps simplifies the overall process, but it complicates error 
determinations, since errors must be propagated through all steps. While there are well-known statistical methods for 
computing confidence ranges in such cases, the many parameters involved (at the very least, one temperature and one 
density for each X-ray spectrum), make the process cumbersome and expensive. Error propagation is more flexible 
when it is possible to marginalize over fitting parameters directly to compute a confidence range for the mass. This is 
considerably easier when the mass can be determined directly in terms of parameters used to fit to the X-ray data, as 
in the approach used here. 

2.1. A model-independent approach 

A simple model can be used to approximate any spherical cluster and its atmosphere well. The cluster volume 
is divided into a series of spherical shells centered on the cluster and, in addition to the assumptions of spherical 
symmetry and hydrostatic equilibrium, it is assumed that 

1) the gravitating matter density in each spherical shell is constant and 

2) the gas in each spherical shell is isothermal. 

Assuming that the gravitating matter density is constant in shells provides, essentially, the simplest continuous mass 
distribution that is nonsingular at r = 0, which can approximate any real mass distribution well with a large number 
of shells. Similarly assuming that the shells are isothermal results in, perhaps, the simplest gas model that can 
approximate any real, spherically symmetric, gas temperature distribution well with a large enough number of shells. 
Projected onto the sky, the spectrum for any region of the model clu ster is a sum o f optically thin, thermal models, a 
form that can be expressed conveniently as an XSPEC mixing model (jXrnaud 19961). 

Although the term is inappropriate in the strict mathematical sense, since this model can approximate arbitrary 
mass and temperature distributions, we call the method "model-independent." The model used here is one instance 
of a broad class of (model-independent and model-dependent) approximations that could be used for the gravitating 
matter distribution (alternative models would differ in equations [2] and [3] below). 

Denoting radii measured from the center of the cluster by r and the shell boundaries by r\,r2, ■ ■ ■ , r n+ i in increasing 
order, the i th shell is the region with rv, < r < fj+i. If the gravitating matter density within the i th shell is pi, then 
the gravitating mass distribution within the shell is given by 

4-7T 

M(r)=M l + — Pl (r 3 -rf), (2) 



where Mj = M(tv,) is the total gravitating mass within rv,. The gravitational potential within the i th shell is then given 
by 

( ' }Ti + ^Gp i (r-r i )(r + 2r i ) 



<r>(r)-<r» (/•,)= / ^^±dr = 



n 3 



^ 0) 



Since the gas in each shell is assumed to be isothermal and hydrostatic, within the i th shell the electron density can 
be expressed as 

^m H {$(r) - $(r 4 )}" 



n e (r) = n e ,jexp 



for rv, < r < n+i, (4) 
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Figure 1. Shells and annuli used for the mass model. This schematic shows a cross-section through a cluster center that contains our 
line-of-sight. Dashed lines show boundaries of cylindrical regions that project onto the annuli on the sky from which X-ray spectra are 
extracted. Full lines show boundaries of the corresponding spherical shells. 

where n e ,i is the electron density at the inner edge of the i th shell, T; is the gas temperature in that shell, $(r) is the 
gravitational potential, k is Boltzmann's constant and fimn is the mean mass per particle in the gas. For the gas to 
be hydrostatic, the pressure must be continuous from shell to shell, but, as the temperature is generally discontinuous 
at shell boundaries, so is the gas density. Pressure continuity between shells i and i + 1 requires 



; exp 



kTi 



Ti. 



(5) 



which determines n e ,t+i in terms of n c ,i, the temperatures and the potential. 

Together, these results show that when the gravitating matter density and gas temperature are specified for each 
shell, under the model assumptions, the gas density is determined throughout the model atmosphere, up to a single 
scale factor that may be regarded as the electron density, n e ,i, at the inner edge of the innermost shell (usually r = 0). 

When observing a spherical cluster, the normal practice is to divide the sky into concentric annuli centered on the 
cluster and extract a spectrum from each annulus. An annulus is the projection onto the sky of a cylindrical region 
in the cluster. The schematic in Fig. [1] shows sections through some cylinders (bounded by dashed lines) in a plane 
that contains our line-of-sight through the cluster center, i.e., the common axis of the cylinders. We now define the 
spherical shells discussed above to have radial bounds matching those of the cylinders (full lines in Fig. [T]). Denoting 
cylindrical radii, measured from the line-of-sight through the center of the cluster, by zu, the j th cylindrical region is 
then defined by Wj < w < zuj+i and, by construction, = vui, which simplifies some of the following results. The 
same geometric arrangement is used for deprojections. 

To complete the spectral model, we need to calculate the emission measure of the gas from the i th shell that appears 
in the j th annulus, i.e., the emission measure of the gas in the intersection between the i th shell and the j th cylinder 
(hatched region in Fig. [T]), 



EM hJ = 



n c nn dV = Air 



Intersection 



i-i+i 



n e (r)n H (r) (yV 2 - r 2 - yV 2 - r| +1 ) 



r dr, 



(6) 



where each square root should be taken as zero if its argument is negative. These integrations are carried out 
numerically. Note that for i < j, since the i spherical shell does not intersect the j th cylinder, EM^j = 0. The model 
spectrum for the j annulus is given as a sum over the shells of the thermal model for each spherical shell, weighted 
by its emission measure in the j th cylinder, i.e., the spectrum is 



F^E) =J2^idf(E,T i ,A i ) ) 



(7) 



„-th 



where f{E,T i ,A i ) is the thermal spectrum for the i spherical shell (for unit emission measure). The parameter, A^, 
is a placeholder for all thermal model parameters, including abundances, apart from the temperature. 

High quality spectra are required to obtain good results when fitting this model. Since the demands of the model 
are similar to those of deprojections, it is expected that the minimum required source count per spectrum is several 
thousand. There is a clear trade off between reducing shell width to gain spatial resolution and increasing it to improve 
spectral signal-to- noise (see section |4~2|) . 



2.2. Implementation 

The model outlined above has been implemented as an XSPEC mixing model called CLMASS. This mixing model 
has two parameters per spherical shell, the gas temperature and the gravitating matter density. The temperature for 
each shell must be linked (equated) to the temperature parameter of a thermal model for that shell. As noted, there 
is a single overall gas density scale that is free in the model. For the mixing model, it is convenient to make the 
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Table 1 

The cluster sample 



Cluster 


z 


^soo 1 


^500 1 


Chandra ObsID 






(kpc) 


(kpc) 




A907 


0.1603 


501 


1095 


3185, 3205, 535 


A1413 


0.1429 


559 


1300 


1661, 5002, 5003 


A1991 


0.0592 


341 


734 


3193 


A2029 


0.0779 


642 


1359 


891, 4977, 6101 


A2390 


0.2302 


561 


1414 


4193 


A1835 


0.2520 


673 


1475 


6880, 6881, 7370 


A1650 


0.0845 


515 


1128 


5822, 5823, 6356, 6357, 6358, 7242 


A3112 


0.0761 


459 


1025 


2216, 2516, 6972, 7323, 7324 


A2107 


0.0418 


416 


919 2 


4960 



a Overdensity radii from Vikhlinin et al. (2006). 

h Chandra data for A2107 do not extend to rsoo, so that M500 cannot be determined for A2107. 

XSPEC norm for the innermost shell the free parameter instead of the central gas density, n e ,i. The XSPEC norm is the 
emission measure multiplied by a scale factor that is the same for each shell and depends on the distance to a cluster. 
In using the CLMASS model, the XSPEC norm for every outer shell is linked (equated) to that of the central shell and 
the X-ray flux for the i th annulus (equation [7]) is given by 

n 

S j (E)=KY,W i jf(E,T i ,A i ), (8) 

where K is the single free norm (for the innermost shell) and there are n spherical shells. The weights, Wjj = 
EMjj/EM^x, are completely determined by the parameters of the mixing model. 

Frequently, a cluster extends beyond the field of view of an X-ray detector, so that there is X-ray emission in the field 
of view originating from gas lying beyond the outer edge of the outermost spherical shell (the n th shell) of the CLMASS 
model. There is no entirely satisfactory way to deal with such emission, but we often want to use such nonideal data 
for mass determinations. Using a local X-ray background measured immediately outside the n th annulus generally 
overestimates the "background" gas emission in all inner shells. Based on empirical fits to gas density profiles, the 
CLMASS model allows the alternative of assuming that the gas density profile beyond the inner edge of the n th shell 
follows a beta model, 

n e (r) =n c , (l + r 2 /a 2 r 3 P/ 2 , (9) 

where f3 and a are constant parameters. The gas density at the inner edge of the n th shell is still determined by pressure 
continuity (equation [5]) , but, instead of using the gas density profile from equation the emission measure for the 
n th shell in each annulus is determined by integrating the beta model to infinity along our line-of-sight (giving results 
in terms of incomplete beta functions). Together with a switch that determines whether or not this feature is used, 
the beta model adds three fixed parameters to the CLMASS model. This feature is at odds with the model-independent 
approach of CLMASS otherwise. Ideally, it should be avoided by using X-ray data that cover the whole of a cluster. 

For the sake of brevity, the new mass model will be referred to as CLMASS in the remaining discussion. Much of that 
discussion would apply to any implementation of the model, although there are some comments on limitations that 
are specific to the XSPEC implementation, as noted. 

3. MASSES FOR A SAMPLE OF GALAXY CLUSTERS 
3.1. Sample and data preparation 

The CLM ASS model has been app lied to Chandra data for a subsample of the clusters analyzed by I Vikhlinin et al.l 
(200<| and iVikhlinin et al.l (p00 9a). Target clusters were chosen to sample a broad range of cluster temperatures, 
hence masses. Clusters analyzed here are lis ted with some o f their basic properties in Tabled] 

To simplify comparison with the results of IVikhlinin et all our analysis parallels theirs closely, using the same X-ray 
spectr a as far as possible with the CLMASS model. Details of the spectral reductions can be found in IVikhlinin et al.l 
(2005). The spectra were extracted from annuli with radial bounds having voj + \ = 1.5n7j, for j > 1. Since the 
gravitating matter density is constant in the corresponding spherical shells, the matter density profile is sampled 
rather coarsely in t he resulting mod el. To make our masses directly comparable, we have also adopted the values of 
r 2500 an d rsoo from IvTkhhnin et al.l (Table [TJ. Note that this eliminates a source of uncertainty in M2500 and M 500 , 
r educing confiden ce ranges compared to ab initio determinations of these masses. 

IVikhlinin et al.l combined spectra from the different front illuminated (FI) chips into single spectra. However, because 
they have significantly different responses, FI spectra were kept separate from those extracted from the back illuminated 
(BI) chips. This means that the number of spectra per annulus for a cluster can vary in the transition regions between 
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Figure 2. Mass profiles for Abell 907. The black line shows the cumulative mass distribution for the best fitting OLMASS model when no 
constraints are applied to the gravitating matter densities. Crosses mark the shell boundaries. Large variations in slope from shell-to-shell 
reflect large variations in the gravitating matter density. Left panel: The cyan line shows the mass profile at the upper 90% confidence 
limit for the mass within rsoo, i.e., M500, and the green line shows the corresponding mass profile at the lower 90% confidence limit. Right 
panel: Same for A^soo- 

FI and BI chips. The version of CLMASS used here inherited from PROJCT a requirement that the same number of 
spectra be used for every data group. Particularly for observations with ACIS-S at the aim point, this means that we 
were unable to use every available spectrum for some annuli. In cases where a choice had t o be made, we o pted for 
the spectrum with the greater source count. As a result, we used slightly fewer spectra than [Vikhlinin et al.1 for some 
mass determinations. This shortcoming of CLMASS will be remedied in a future version^ 

A limitation of XSPEC mixing models is that they must be applied after all other model components. This is 
inconvenient for several reasons. Photoelectric absorption due to foreground gas should be applied to the spectrum 
emerging from each annulus, after it is emitted by a cluster. In practice, the foreground absorption usually does not 
vary significantly across a cluster, so that it is satisfactory to apply the same absorption model to each thermal model 
before the mixing model. However, there are clusters where the absorbing column varies significantly, which cannot 
be modeled well by the CLMASS model in xspec. 

Another inconvenience arises in background corrections. In the standard approach, background spectra are con- 
structed from source free exposures, processed and scaled to match the observations. Variations of the particle 
background and the soft X-ray background with time and position can leave residual backgr ounds requiring furthe r 
correction. These are often addressed by model components added to the annular spectra (| Vikhlinin et alJ T2005) . 
but this approach cannot be used in conjunction with an xspec mixing model0 Here, the residual backgrounds are 
corrected using a specially prepared CORRFILE for each annular spectrum. Depending on need, one or two soft thermal 
components were fitted to the spectrum from the outermost annulus for each cluster (to account for a residual soft 
background and any emission from the North Polar Spur). These components were scale d by area and added to the 
existing CORRFILE for each annulus (already employed to correct for the read out artifact: I Vikhlinin et al.ll2005D . 

Gratifyingly, inclusion of the beta model feature described in section 12721 had little effect on the quality of the model 
fits and on the masses. This is to be expected, since the clusters are barely detected in the largest annuli, so there 
should be little need to account for cluster emission from regions beyond the outermost shells. 

3.2. Cluster Masses 

The CLMASS mixing model was combined with an absorbed thermal model, W ABS*MEKAL, to fi t the sample clusters. 
Column densities for photoelectric absorption were set to the values used by [Vikhlinin et al.l (|2005| ) . As discussed 
in section 12. 2[ many parameters must be linked to use the CLMASS model. Since the number of shells provided by 
the model is fixed, usually there are unused parameters that must also be frozen. Parameters were given reasonable 
initial values to forestall numerical problems ( section 12.21) . With ~ 10 annuli per cluster, in order to avoid error 
prone manipulations of numerous parameters, these steps were automated in xspec Tel scripts. Abundances were 
generally free in the fits, but where they were poorly constrained or took on unreasonable values, abundances for 
groups of adjacent shells were tied together to produce physically reasonable results. Since they have little impact on 
the masses, we do not consider abundances further. All free parameters and links used for the fits were maintained 
when determining confidence ranges for the masses. 

5 A new version of the XSPEC model addressing this and some other issues has been completed. The new version also makes it easier 
to use alternative forms for the gravitational potential and a model using the Navarro, Frenk &; White potential has been added as a 
demonstration. Verions of these models are also available for CIAO SHERPA. 

6 In XSPEC12, such components can be added by defining additional models. 
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Figure 3. Mass profiles for Abell 907 with the gravitating matter density constrained to be monotonic. Details as in Fig. [2] 
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Figure 4. Temperature profiles for Abell 907. Bla ck triangles show tem peratures for the best fitting "monotonic" CLMASS model. Blue 
crosses show the (projected) temperature profile of Vikhlinin ct al. (2005). Temperature errors are 68% confidence ranges for individual 
shells. Radial error bars show extents of the annuli. CLMASS temperatures are deprojected, which generally makes them noisier than 
projected temperatures. 



Initial attempts to fit the data showed that the CLMASS model does not generally constrain the gravitating matter 
density tightly for individual shells. This is because the gravitational potential is governed by the accumulated mass 
within any radius and variations of the gas density depend on changes in the potential. Since there are two integrations 
between the mass density and changes in the potential, the gas density, hence the spectra, are only affected at second 
order by local changes in the gravitating matter density (equation [3]) . Splitting the gravitating matter from one shell 
into the two adjacent shells, in suitable proportions, has little impact on the fit. In the presence of noise, this freedom 
to exchange mass from shell to shell results in large shell to shell variations of gravitating mass density, producing 
noisy best fitting mass profiles like that for Abell 907 in Fig. [21 

The freedom of gravitating mass to shift between adjacent shells can be limited to some extent by constraining 
the gravitating matter density to be a monotonically decreasing function of the radius. This constraint improves 
model behaviour with a relatively minor impact on the goodness of fit (Table [2]). In practice, the constraint requires 
gravitating mass densities to be linked for groups of adjacent shells where the unconstrained density would otherwise 
increase with radius. Mass profiles for Abell 907, with the gravitating matter density forced to be monotonic, are shown 
in Fig. [3l Evidently, the monotonic constraint makes the matter density profile considerably smoother. Although this 
condition is not required of the gravitating matter density distribution, it is simple, physically reasonably and relatively 
weak, while it reduces the noise in the mass profiles significantly. As a consequence, we prefer results for "monotonic" 
mass profiles in the following, although some results for unconstrained mass profiles are also given for comparison. 

Once the temperatures have been determined, gas densities are governed by the distribution of gravitating matter. 
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Table 2 

Cluster masses 



Unconstrained 



Monotonic 

Aftoo M2500 
(10 14 Mq) (10 14 M ) 



Cluster 



X 2 /dof M500 

(10 14 Mq) 



X 2 /dof 



A907 

A2390 

A1835 

A1650 

A3112 

A2029 

A1991 

A1413 

A2107 



692.4/626 3.2; 

840.9/626 23.0 

489.4/290 7.1 j 

669.1/301 5.1; 

545.8/314 3.0; 

2988.0/664 6.6; 




693.7/630 
846.2/632 
491.7/293 
673.2/306 
546.7/317 
2993.5/669 
563.7/326 
450.3/321 
428.1/299 



4.7] 
14.9 

Ml 
5.8^ 
2.8 H 

r.ij 

5.9" 1 
9.lj 




449.0/313 



7.9 



,+ 10.5 
-3.9 



Loosely speaking, when fitting the CLMASS model, temperatures are determined by spectral shapes, so that the 
gravitating matter densities are determined by the radial distribution of the X-ray emission (i.e., emission measure). 
However, this is only a rule of thumb. Especially in the outskirts of clusters where the X-ray emission is not well detected 
above the background, fitted temperatures may be poorly constrained by the spectra, so that the temperatures can 
also respond to variations in X-ray brightness. This may result in large temperature excursions in response to poorly 
modeled fluctuations in X-ray brightness. Also, since the temperatures apply to gas in the spherical shells, they are 
deprojected temperatures. Like other deprojected temperatures, the spectral models are coupled from shell-to-shell, 
making the temperatures etc., prone to fluctuate from shell-to-shell. This effect can be seen in Fig. [H which compare s 
the best fitting CLMASS temperature profile for Abell 907 to the projected temperature profile of Vikhlini n et al.l ((2005). 

Cluster masses obtained from the CLMASS model are given in Table [2j Consistent with the model assumptions, 
gravitating matter densities are treated as constant within each shell, so that mass profile within a shell is as given by 
equation ([2]) . Best fitting chi squareds for the unconstrained fits are given in column 2 and the corresponding values of 
M500 are in column 3. The best fitting chi squareds for the monotonic fits are in column 4, with corresponding values 
of Af5oo and -M2500 in columns 5 and 6, respectively. Values used for rsoo and r25oo are given in Table [TJ There are no 
results for M^q for Abell 2107 because the Chandra data for it do not reach r^. No unconstrained result is given for 
Abell 1991, the coolest and faintest cluster in the sample, because the unconstrained CLMASS model places no useful 
90% upper limit on M500 for it. Abell 1991 is not detected at high significance in annuli beyond rsoo, so that the 
emission measures of the corresponding shells can be made arbitrarily small in the models. This permits an arbitrarily 
large gravitating mass in the shell immediately inside r$oo- Making the gravitating matter density monotonic prohibits 
the density in the outermost shell within rsoo from exceeding that in the next inner shell, so that the mass of Abell 1991 
can no longer be arbitrarily large. Thus we do get a confidence range for -M500 using the constrained CLMASS model. 

Based on chi squared, all of the fits are formally unacceptable. However, all clust ers show some departures from 
spherical symmetry and many also show evidence of ongoing minor merger activity (e.g. lMarkevitch & Vikhlininll2007h . 
so that we should expect modest departures from hydrostatic equilibrium too. Considered in this light, it is gratifying 
that the coarsely resolved CLMASS model used here gives reduced chi squareds of less than 2 in most cases. The clear 
exception is Abell 2029, with reduced chi squareds of ~ 4.5. Abell 2029 is the only cluster in the sample for which there 
are both ACIS-I and ACIS-S observations, giving two complete sets of spectra. Fitting these separately gives reduced 
chi squareds of 1.4 for ACIS-I and 2.6 for ACIS-S, more consistent with the other results. The data for ACIS-I are a 
lot shallower than those for ACIS-S (exposures of ~ 10 ksec and 100 ksec, resp.), which is why the reduced chi squared 
appears better for ACIS-I. However, the chi squared for the joint fit is more than a factor of 2 larger than the sum of 
the chi squareds for the separate fits, pointing to incompatibility between the two sets of spectra. Most likely, this is 
due to departures from spherical symmetry in Abell 2029. Especially on larger scales for the ACIS-S detectors, the 
field of view is restricted, so that spectra are not sampled uniformly from every annulus. Since the regions sampled by 
ACIS-I and ACIS-S differ, anisotropies can lead to significant differences in the two sets of spectra. This effect is only 
evident when there is more than one set of spectra for a cluster. It highlights the imperfections of our assumptions. 



To find confidence ranges for the mass, the mass is varied about its best fitting value and, for each trial mass, a new 
best fitting model is determined. The mass range for which the best fit statistic remains below an appropriate limit 
then defines the confidence range for the mass. One virtue of the CLMASS model is that calculation of mass confidence 
ranges can be undertaken as part of the fitting process, making it possible to marginalize over all free parameters of 
the fit. 

Fixed mass constraints are straightforward to implement. The total gravitating mass within the radius r is simply 
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Figure 5. Cluster masses, M500, from IVikhlinin et~ al, (2005) plotted against those from the OLMASS model. The unconstrained best fits 
are shown in blue, with their 90% confidence ranges in green. With the monotonic condition ap plied, the best fits are shown in red and 
their 90% confidence ranges in black. 90% confidence ranges are also shown for the masses from|Vikhlinin ct al. (2005). The dotted line 
shows equality. No unconstrained result was obtained for Abell 1991. 

where 

Vi{< r) = — <r 3 - rf, for r t < r < r l+1 (11) 
1 0, otherwise, 

is the volume of the i th spherical shell lying inside r. For a fixed mass, M(r), provided that Vk(< r) ^ 0, this is 
equivalent to 

_ M{r) v V t (<r) 
V k (<r) V k (<r) 

which may be used as a constraint on pk- If Pk bumps up against a further constraint (e.g., it goes to zero, or 
it becomes non-monotonic when the monotonic condition is applied), that constraint must take precedence. The 
fixed mass constraint can then be applied to another shell. In cases where densities of adjacent shells have to be 
tied to maintain the monotonic constraint, the fixed mass constraint can be applied to a group of shells with the 
same gravitating mass density (in equation II 2\ Vk{< r) is replaced by the total volume inside r for the group of 
tied shells and the sum must exclude all shells of that group). In general, managing constraints that may conflict 
with one another is complex. In an attempt to minimize such conflicts, the fixed mass constraint was applied to 
the shell with the largest volume, Vk(< r), since that choice makes p k least sensitive to pi (under the constraint, 
dpk = -T, i7 tk d Pi V i(< r )/ y k(< r)). 

Much the same issues apply for determining confidence ranges for the mass at a fixed overdensity, rather than at a 
fixed radius. Specifying the overdensity determines the relationship between r and M{r) (i.e., M(r) — 47rpr 3 /3, where 
the overdensity, p, is fixed), so that r and Vi(< r) vary as M varies. However, since r remains fixed for a fixed mass, 
the remainder of the discussion above is unaltered. 



4. DISCUSSION 
4.1. Comparison to conventional mass determinations 

Fig. [5] compares values of M500 determined from the CLMASS model to those obtained by IVikhlinin et al.l ((2006) 
using essentially the same spectra. The figure shows 90% confidence ranges (A% 2 = 2.706) for both the case where 
the gravitating matter density is unconstrained (best fit blue, confidence range green) and when it is constrained to be 
monotonically decreasing (best fit red, confidence range black). Plainly, the CLMASS model does not place tight limits 
on values of M500 for the sample clusters. For the unconstrained fits, the 90% confidence intervals typically extend 
over a factor of ~ 3.5 in mass. 

The situation improves when the monotonic constraint is applied. Apart from Abell 1991, the average length of the 
90% confidence ranges is a factor o f a 2.3 in mass ( averag ed in log space), compared to ~ 1.34 (i.e., SM/M ~ 30%) 
for the model-d ependent results of IVikhlinin et all ((2005). The formal agreement between the CLMASS results and 
IVikhlinin et al.l is good, apart from Abell 1991. The CLMASS value of M500 is high for such a cool cluster, stemming 
from a temperature for the outermost shell inside rsoo that is high (kT ~ 8 ± 4 keV, la, compared to a projected 
temperature of ~ 2 keV). With a sample of 8 clusters, it is reasonable to expect disagreement at the 90% level in 
one case. However, the poor result for Abell 1991 mainly reflects inadequate data at rsoo and beyond. This result 
highlights the more general issue for the CLMASS model, that the total mass is only well-constrained if we have good 
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Figure 6. Cluster masses, 

A%00, from Vikhlinin ct al. (2005) versus those from the OL MASS model. Best fittin g results from CLMASS with 
the monotonic condition and their 90% confidence ranges are plotted in x and results from Vikhlinin ct al. (2005) with their 90% confidence 
intervals are plotted in y. The dotted line shows equality and the full line shows the best fitting straight line. 

X-ray data for annuli inside and outside the radius of interest. The lack of better data from beyond r5oo is the primary 
cause of the large 90% confidence intervals for M 50 o. 

Masses are better determined at smaller radii, where we have good data from beyond the radius of interest. For 
the sample clusters, r 2 soo i s roughly half o f .ggpo (Table [D. Values of M2500 from monotonic CLMASS model fits are 
compared to the results of Vikhlinin ct al. (2005) in Fig. [6] The results for M2500 are markedly more accurate than 
those for M 500 , clearly demonstrating the benefit of having high quality spectra that bracket the radius of interest. 
Reflecting the greater freedom of the model, 90% confidence ranges from CLMASS remain larger than those for the 
model-dependent masses, with a mean spread for the sample of 45 %, compared to 25% for the model-dependent 
masses. Individual values of M2500 also agree well with the results of IVikhlinin et all (I2005D . 

There is a modest off s et be tween the CLMASS results overall and those of IVikhlinin et al.l Using the bisector method 
of lAkritas fc Bershadvl (|1996l ) to fit the relationship between the masses 
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gives a = —0. 86 an d b = 1.101, which is plotted in Fig. |B1 Monte Carlo simulations using the masses from 
IVikhlinin et al.l (|2005l ) and the fitted errors from both methods show that b is within one standard deviation of 
unity However, at the sample mean (in log space), the fitted relationship gives masses that are 12% low (i.e., the 
CLMASS results are high), which occurs b y chance about 2 % of the time. Thus, there is a moderately significant offset 
between the values of M2500 obtained by IVikhlinin et akl and those obtained here. 
As noted above, the annuli used here are wider than ideal, so that the mass profile used by the CLMASS model is 

coarsely sampled. The gravitational potential, A$ = J b GM(r)/r 2 dr, is more sensitive to mass in the inner part of 
a shell than the outer part, but, relative to a matter (density distribution that decreases with radius, the constant 
density approximation used by CLMASS puts more mass in the outer part. Thus, slightly more mass is required in the 
coarsely sampled CLMASS model to obtain the same change in gravitational potential. Since the distribution of gas 
density (i.e., X-ray emission) is determined by the gravitational potential, that is the physical quantity which is fixed 
by fitting the X-ray data. As a result, the coarsely sampled CLMASS model finds masses that are biased slightly high. 
An estimate of this effect is derived in the Appendix and the bias is found to be no more than ~ 3%. Allowing for the 
bias expected for r\ = 2 (see Appendix), the likelihood of the observed 12% bias at the mean mass is increased from 
about 2% to ~ 7%. We note that the bias due to coarse sampling is quadratic in (r^+i/r^ — 1), so that it decreases 
rapidly with finer sampling. The significance of the remaining ~ 9% mass bias is marginal. The issue of possible bias 
in the CLMASS model will need to be addressed by other means. 



4.2. Strengths and Weaknesses 

The CLMASS model can approximate arbitrary distributions of temperature and gravitating matter, but this freedom 
means that it fails to constrain cluster masses as tightly as other methods of X-ray mass determination employing more 
restrictive assumptions. Gravitating mass can also move relatively freely between adjacent shells in the CLMASS model 
with little impact on the fit. Its performance in this respect can be improved significantly by requiring the gravitating 
matter density to be a monotonically decreasing function of the radius. This results in smoother mass profiles and 
reduces confidence intervals for the mass. The effect of adding a relatively weak assumption to the CLMASS model 
illustrates the trade offs between the (apparent) precision of mass determinations and the assumptions they rely on. 
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Using some of the best X-ray data available and insisting that the mass distribution be monotonic, the CLMASS 
model still fails to limit M^qq to within a factor ~ 2 at the 90% confidence level (Fig. [5]). However, it gives well- 
constrained results for Af25oo- The poor results for M500 are a consequence of the generality of the model and the 
relatively poor quality of the X-ray data near and beyond rsoo ■ This situation will be slow to improve, because X-ray 
emission is proportional to the square of the gas density, making it decline rapidly at large radii. It suggests that 
other X-ray mass determinations rely heavily on their model-dependent assumptions to obtain accurate masses at 
large radii. While the issue is most acute for X-ray data, mo del- dependent assumptions have the same effect on all 
mass determinations. A mo del- dependent prior, such as the assumption that a cluster potential has the NFW form, 
reduces confidence intervals for its mass. Unless a prior is accurate, confidence ranges for the mass that rely on it will 
be too optimistic. Discrepancies between mass determinations may be caused by inappropriate priors as much as by 
failures of the physical assumptions that underlie mass determinations. 

Except for cl usters that are obv iously dynamically active, X-ray mass determi nations generally agree well with dy- 
nami cal masses (Rines et al. 2003|) and with masses determined from weak lensing (jMahdavi et al.| [2008; Me neghetti et al.l 
I2009T) . with a modest bias due to departures from hydrostatic equilibrium. Note that the sample of iMahdavi et al.l 
( 2008) includes a number of dynamically active cluster s, so that their mean bi as of ~ 20% in M500 overestimates the 
bias for relaxed clusters, such as those used here and bv I Vikhlinin et aD (j2009b[ ). Since mo del- dependent X-ray masses 
agree as well as expected with other mass determinations, the confidence ranges for M500 obtained with the CLMASS 
model appear unnecessarily poor. On the other hand, using the excellent X-ray data bracketing r25oo, the CLMASS 
confidence ranges for M2500 are not a lot more conservative than those obtained by careful modeling of the X-ray data 
(|Vikhlinin et al.ll2005l) . Until there are better X-ray telescopes, use of the CLMASS model is probably better restricted 
to the brighter regions of clusters. The results here indicate that mo del- dependent X-ray determinations of M 500 rely, 
in effect, on extrapolating to rsoo models that are largely constrained by data at smaller radii. In view of this, using 
M2500 obtained by fitting the CLMASS model as a proxy for the total mass, also an extrapolation in effect, may serve 
as well. 

As noted (section l3.ip . the annuli used here have zuj+i/zuj = 1.5, making the gravitating matter density coarsely 
sampled for the CLMASS model. While spectra with many photons are needed to obtain accurate temperatures (more 
so because CLMASS effectively deprojects the temperatures), binning the data coarsely discards information about 
the radial distribution of surface brightness that could help to constrain the distribution of gravitating mass. The 
performance of the CLMASS model may well be improved by extracting spectra from narrower annuli. It may even 
be advantageous to bin the data so finely that individual spectra alone have insufficient counts to constrain the 
temperatures well. Spectral parameters for groups of adjacent shells could be tied, in effect treating them as a single 
shell for the spectral fits, while exploiting the information from the individual spectra on the distribution of brightness. 
This and related issues, such as the best choice of fit statistic, remain to be explored. 

There is little constraint on the thermal models that can be used for the spherical shells with CLMASS. In particular, 
it would be feasible to use multiphase models (e.g. iRussell et aLll2008t ). provided that an appropriately weighted mean 
temperature is used for the CLMASS model temperature for each shell. If the multiphase gas in a shell is well-mixed 
and in local pressure equilibrium at pressure p, the appropriate temperature, T, can be determined from the mean gas 
density, p g , by kT = iimnp/~p g . In terms of the properties of the separate phases, this yields the weighted temperature 

T = (j2 q K q T^j I (j2 q KqTq^j , where K q is the xspec norm of the q th phase in the shell, T q is its temperature and 
the sum runs over all gas phases in the shell. 

5. CONCLUSIONS 

A new approach for making X-ray mass determinations of galaxies, groups and clusters of galaxies was introduced. 
The new method relies on the usual assumptions of hydrostatic equilibrium and spherical symmetry. Cluster gas is 
also assumed to be isothermal and the gravitating matter density constant in spherical shells centered on a cluster. 
With sufficient shells, this distribution provides a good approximation to any spherically symmetric distribution of 
temperature and gravitating matter, so the approach is largely model- independent. It does not require the additional, 
model-dependent assumptions that are essential to most conventional X-ray mass determinations. 

The method has been implemented as an XSPEC mixing model called CLMASS. With no further constraints, the 
model only determines the gravitating mass in individual shells poorly. This behavior is improved significantly by 
constraining the gravitating matter density to be a monotonically decreasing function of the radius. 

The CLMASS model w as used to determine values for M500 and M2500 for nine relaxed clusters from the sample 
of I Vikhlinin et all (|2005| ). as far as possible, using the same spectra. While the results for M 50 o agree with those of 
iVikhlinin et al.l . the CLMASS confidence ranges are much larger than those for the more conventional method. This 
reflects the lack of assumptions about the form of the gravitating matter distribution in CLMASS and a lack of good 
X-ray data from regions beyond ^00- At smaller radii that are bracketed by good X-ray data, the CLMASS model 
pe rforms much be t ter. F or M2500, CLMASS gives confidence ranges that are only moderately more relaxed than those 
of IVikhlinin et al.l (|2005f >. Individual masses agree well between the two methods, but there is a small offset in the 
relationship between the two determinations of M2500 at the mean mass. This is partly a consequence of the coarse 
radial sampling of the spectra used here, but there is a marginally significant offset beyond that. 

In regions where we have good X-ray data, the new method can determine masses with comparable accuracy to 
conventional methods, without relying on model-dependent assumptions. Until we have better X-ray data for the 
faint outer regions of clusters, it will be less useful for determining total cluster masses, other than by using M2500, 
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for example, as a proxy for the virial mass. In contrast to most other approaches, the new method allows confidence 
ranges for the mass to be determined directly from fits to X-ray spectra. The method can be used with a wide range 
of thermal models for the X-ray emission from individual shells within a cluster, including multiphase models. 
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APPENDIX 

COARSE SAMPLING MASS BIAS 

For the purpose of estimating the mass bias dues to coarse sampling of the distribution of gravitating matter, we 
consider a model cluster in which the density distribution is a power law, pt{r) oc r~ r ', with rj < 3. Locally, this 
serves as a reasonable approximation for more realistic mass distributions. As argued in section 231 the gravitational 
potential determines the gas distribution, so we assume that the potential is most constrained by the data. Thus we 
need to compare potentials for the "true" cluster and the CLMASS approximation to it. Corresponding to the power 
law density distribution, the mass within radius r is 

M t (r) = M t 4r/r t f-\ (Al) 

where Mt,i is the mass within rj. For this mass distribution, the change in gravitational potential from n to is 

J r . r 2 n 2-T] 

where 1 + A = r i+ i/ri. Note that for 77 — 2 , the expression [(1 + A) 2 ~ r ' — l]/(2 — rj) — >• ln(l + A) is well behaved. 

Now consider the CLMASS model approximation for this cluster. As for the sample data, we assume that the shell 
boundaries follow a geometric progression, r i+1 = (1 + A)r*j, for all i, but we will ignore the cutoff at i = 1, assuming 
that the geometric progression continues to i — — 00 to simplify the sums (an accurate approximation). Matching the 
power law for the true density distribution, we also assume that the mass densities for the CLMASS model follow the 
geometric progression 

Pi+i/Pi = (n+i/ri)"" = (1 + A)-". (A3) 
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This condition is not strictly required, but, like the power law for the "true" density, it should be a reasonable local 
approximation. The form (IA3|) determines the ratio AMj/M^ found below (IA5|) . making it critical to our estimate of 
the mass bias. For the CLMASS model, the mass in the i th shell is then 



AM, = yft(rf +1 - r?) = y ft rf[(l + A) 3 - 1]. 



(A4) 



Under our assumptions, we have AM& = AM,*(1 + A)^ ^ , so that the shell masses form a geometric progression. 
Thus, the total mass within t{ for the CLMASS model can be summed to give the total mass within 7"j, 
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(A5) 



The increase in potential in the i th shell is (equation [3]) 



A$ ?; = 



GM, 2tt 



— G Pi rfA{3 + A) 



and using (IA4[) and (| A5|) to express Gpirf in terms of GM i jr il this becomes 

A(3 + A) (l + A) 3 -"-l 
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(A7) 



If the potential is determined exa ctly by th e data, as argued above, then fitting the CLMASS model should make 
A$j = A$ t i . From equations (|A2[) and (|A7|) . this determines the ratio of the true mass to the mass obtained from 
the CLMASS model as 
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(A8) 



where 6 is the mass bias. The last form is an approximation showing the lowest order deviation from unity. The 
CLMASS model makes the matter density uniform in shells, so there is no bias for a uniform density distribution, i.e., 
for 77 = 0. The bias is also zero for 77 — ^ 3, since then the mass is concentrated at the cluster center and A$ 4 is 
dominated by the potential due to mass inside fj, which is determined correctly in CLMASS. For the cluster sample 
used here, A = r,;+i/ri — 1 = 0.5. The bias for A = 0.5 is overestimated somewhat (~ 60%) by the approximation 
b ~ 77(3 — ?7)A 2 /12. Nevertheless, 6(A = 0.5) is maximized for 77 = 1.5 at b ~ 0.030. Observed values for 77 at r25oo are 
about 2 or a little greater. For 77 = 2, the model described here gives b ~ 0.026. 



